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Abstract. CO (J = 2-1) and SiO(J = 5-4) emission has 
been observed from the molecular envelope around the 
nearby S star tt 1 Gru. The CO line profile differs from 
the usual parabolic shape seen in uniformly-expanding en- 
velopes; it has a Voigt-like profile and two horns. A model 
for line formation in the envelope shows that a tilted, ex- 
panding disk reproduces the observations well. The star 
also has a fast molecular wind, with a projected outflow 
speed of at least 70 and perhaps as high as 90 km s _1 . 
The fast wind is presumably ejected from the poles of the 
disk. These observations show that the complex structure 
seen in many planetary nebulae, including quadrupolar 
structure and fast winds, may largely evolve from struc- 
ture formed while the progenitor star is in the last stages 
of evolution on the AGB. 

Key words: radio lines: stars — stars: AGB and post- 
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1. Introduction 

Mass loss from stars on the Asymptotic Giant Branch 
(AGB) takes place as a cool, dusty, molecular wind at 
rates which are so high that they dominate the evolution 
of the star. The energy in these huge stars is produced 
in their very small cores, leading to the expectation that 
the structure of the star and of its circumstellar outflow 
will be spherically symmetric. However, a plethora of re- 
cent observ ations (e.g. Lucas 1999 Monnier 1999 ; Tuthill 
et al. 1999| Monnier et al. 1999 ) have shown that the out- 
flows are far from smooth, steady, uniform and spherically 
symmetric, and that asymmetries can be seen on a wide 
range of size scales. This paper discusses molecular line 
observations of the envelope of tt 1 Gru, which is known to 
have complex structure and large-scale asymmetry from 



previous CO observations (Sahai 1992, hereafter S92) 



tt Gru is a nearby mass-losing S star, at a distance of 



153 pc (Perryman et al. 1997; van Eck et al. 1998). S92 



ing stellar winds. Further, the envelope has north-south 
kinematic structure, and at least two velocity components; 
a "normal" outflow with an apparent speed of about 11 
km s _1 , and a faster outflow with a speed of > 38 km s — 1 . 

The molecular line observations measure a mean mass 
loss rate of several xlO -7 M Q yr _1 , while the 60/im emis- 
sion from circumstellar dust shows that the envelope is ex- 
tended (Young et al. 1993), with a radius of about 6 x 10 17 
cm and a kinematic age of ~ 10 years, tt 1 Gru is the pri- 
mary of a wide binary system, whose secondary, at a pro- 
jected separation of 2.7", appears to be a solar-mass main 



sequence star (Proust et al. 1981; Ake & Johnston 1992). 
The stars are thus at least 6 x 10 15 cm apart and, since the 
system mass is at least 2 M©, the orbital period is about 
6000 years and the orbital velocity about 2 km s _1 . 

Given the possible presence of a fast molecular wind 
and bipolar structure in the envelope of tt 1 Gru, plus the 
presence of its binary companion, we reobserved the enve- 
lope with improved sensitivity and modeled its line emis- 
sion in an attempt to better define its geometry and kine- 
matics. These observations are described in the next sec- 
tion, and the results are given in Section 3. Section 4 de- 
scribes a model of the structure of the circumstellar enve- 
lope. These results are discussed in Section 5, and Section 
6 contains the conclusions. 



2. Observations 

The properties of tt 1 Gru from the literature, plus those 
derived in the present paper, are summarized in Table 1. 
The distance is from the Hipparcos parallax (Perryman et 
al. 1997 ), the luminosity from van Eck et al. ( 1998 ), and 



found that the CO line profiles have a very unusual shape: 
they are double-horned, and do not have the steep sides 
typical of molecular line emission from uniformly expand- 



the variability type and pe riod from Proust et al. ( 1981 ), 
and Kerschbaum & Hron ( 1994 ) . 

The observations were made with the 10.4m Robert B. 
Leighton telescope of the Caltech Submillimeter Observa- 
tory on Mauna Kea, Hawai'i, on September 4-6 1997. The 
southerly declination of the star means that it is always 
at a low zenith angle as observed from the CSO, and the 
consequent high air mass makes submillimeter observa- 
tions, which have higher angular resolution, impractical. 
The star was observed in the relatively transparent 1.2 
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Table 1. Observed Properties of 7r 1 Gru 



q(1950) 

<5(1950) 

Spectral Type 

Variable Type 

Period 

Distance 

Lboi 

SiO(6-5): 
I 

Tmb 

Vo 
V c 

CO (2-1) (high resolution) 

Projected V 

V c 

Mean M 

CO(2-l) (map): 
Projected V 
V c 
M 

Vo 



22 h 19 m 41.13 s 

-46° 12' 02.4" 

S5.7 

SRb 

150 d 

153 +28, -20 pc 
7450 L Q 



0.75 ± 0.05 K x km s _1 
0.042 ± 0.005 K 
13.1 ± 2.1 km s _1 
-12.6 ± 0.8 km s" 1 



5.1 km s^ 1 

-11.8 km s" 1 

4.5 x 10 -7 M yr" 



6.4 km s" 1 
-11.7 km s _1 
1.2 x 10~ 6 M yr" 
15 - 18 km s _1 



CO(2-l) line with both the 50 MHz (high resolution) and 
500 MHz (large bandwidth) spectrometers. The high reso- 
lution observation was made to measure the detailed shape 
of the strong CO line emission close to the systemic ve- 
locity of the star, while the large bandwidth observation 
was made to measure the weak emission from the star's 
high velocity molecular outflows. During these observa- 
tions, the telescope pointing was checked roughly once per 
hour using the strong line emission. Second, we observed 
the SiO(6-5) line; because this line requires high density 
for excitation, and because much of the Si condenses onto 
dust close to the star, this emission arises from the inner 
regions of the envelope. 

The third observation was a map of the CO(2-l) emis- 
sion from the envelope made on a 100" x 100" square grid 
centered on the star's position and sampled every 10". 
The observation was made using the 'on-the-fly' (OTF) 
technique. This method consists of scanning across the ob- 
ject's position and accumulating spectra as the telescope 
moves. The grid is built up by a set of scans stepped in 
the orthogonal direction. The map of it 1 Gru was made 
by executing several rapid complete OTF maps and sum- 
ming them; this procedure minimizes errors due to point- 
ing drifts. The OTF observations were made with the 500 
MHz spectrometer only. Before each OTF map, the point- 
ing was updated using the CO(2-l) line emission. 



mm window, in the CO(2-l) line at 230.538 GHz and the 
SiO(6-5) line at 260.518 GHz. The telescope half-power 
beamwidth was measured to be 30", and the main beam 
efficiency to be 69%, using observations of Jupiter. The 
zenith atmospheric opacity r Q was ~0.06. 

The observations used a liquid helium cooled SIS junc- 
tion receiver with a double sideband system temperature 
of ~100 K. The temperature scale and the atmospheric 
opacity were measured by comparison with a room tem- 
perature load. The temperature scale was corrected for 
the atmospheric opacity and for the main beam efficiency. 

The spectral lines were observed using two acousto- 
optic spectrographs with bandwidths of 500 MHz and 50 
MHz respectively over 1024 channels. The spectrometer 
frequency, frequency scale and spectral resolution were 
calibrated using an internally generated frequency comb, 
and the velocity scale was corrected to the Local Standard 
of Rest (LSR). The velocity resolutions of the spectrome- 
ters were ~ 2.1 km s _1 and ~ 0.21 km s _1 . 

The observations were made by chopping between the 
star position and an adjacent sky position with the sec- 
ondary mirror, using a chop throw of 120" in azimuth at 
a rate of 1 Hz. Pairs of chopped observations were made 
with the star placed alternately in each beam. The spec- 
tral baselines resulting from this procedure are linear to 
within the r.m.s. noise for both spectrometers. 

We performed three observations of the ir 1 Gru en- 
velope. First, we obtained a sensitive observation of the 



3. Results 

The individual observations of the SiO(6-5) line were ex- 
amined for bad baselines, co-added, and multiplied by the 
main beam efficiency. The resulting total SiO(6-5) line 
profile is shown in Fig. [l] (where it is compared with the 
high resolution CO(2-l) line profile, see below). The ob- 
served properties: the integrated intensity in K x km s _1 ; 
the peak line temperature; the wind outflow speed V ; 
and the central velocity V c , are listed in Table 1. The cen- 
tral velocity and outflow velocity are in good agreement 
with the values found by S92 for the 13 CO(1-0) line. The 
SiO(6-5) line, like the 13 CO(l - 0) line observed by S92, 
is roughly parabolic in shape, although the signal-to-noise 
ratio is not sufficient to definitively rule out a more com- 
plicated shape, such as the presence of horns. The profile 
does, however, show that the outflow velocity of the dense 
gas in the inner envelope is about 13 km s _1 . 

Fig. [ljalso shows the high resolution CO(2-l) line pro- 
file. The shape of this line is quite different from that of the 
SiO(6-5) line profile: it has two horns and a non-parabolic 
Voigt profile shape with a larger velocity extent. The mean 
velocity of the horns is -11.8 km s _1 , in good agreement 
with the central velocity of the SiO(6-5) and 13 CO(1-0) 
lines. The horns arc offset in velocity from the line center 
by ±5.1 km s , about half the outflow velocity given by 
the SiO(6-5) line. 

Fig. U shows the broad-band CO(2-l) line profile. 
High velocity molecular line emission is seen, extending 
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Fig. 1. CO(2-l) and SiO(6-5) line profiles from the cir- 
cumstellar envelope of tt 1 Gru. Dark line: CO(2-l) line 
profile observed at a resolution of 0.2 km s _1 . Light line: 
SiO(6-5) line observed at a resolution of 2 km s _1 . Dotted 
line: least-squares parabolic fit to the SiO(6-5) line pro- 
file. The temperature scale of the SiO(6-5) line has been 
multiplied by 20. 



from about -100 km s _1 to at least +40 km s _1 , and per- 
haps to +80 km s _1 . This velocity range corresponds to 
-88 km s _1 to +52 km s _1 or + 92 km s _1 with respect to 
the central velocity. The total extent of 140 — 180 km s _1 
is twice that seen in the lower signal to noise ratio obser- 
vations of S92, and is almost certainly a lower limit to the 
true velocity extent of the fast wind. 

Fig. H shows the CO(2-l) line profiles in a region ±50" 
with respect to the star's position. The emission is only 
slightly resolved by the 30" beam, but kinematic structure 
is clearly seen. This structure was examined by construct- 
ing a set of 11 x 11 pixel (10" spacing) channel maps, i.e. 
maps of the emission at each velocity. The position of the 
emission peak in each map was measured using the routine 
IMFIT from the National Radio Astronomy Observatory's 
AIPS data reduction package to fit a two-dimensional el- 
liptical gaussian model. The results are shown in Fig. [|, 
which shows the declination and right ascension of the 
emission peak at each velocity. Near the systemic veloc- 
ity, where the emission is strong, the data are plotted for 
every channel. At higher relative velocities, the data are 
averaged over 10 - 15 velocity channels. Also indicated in 
Figure 4 are the velocities of the emission horns deter- 
mined from the central profile (Fig. [l]) and the velocity 
range (V c ± V Q ) of the SiO emission (Fig. |l|). 



Fig. 2. CO(2-l) line profile of tt 1 Gru (heavy line). Light 
line: scale increased by a factor of 10. The expanded profile 
has also been binned to a resolution of 4 km s . Insert: 
the average CO (2-1) line profile from mapping observa- 
tions over ±30" about the star's position (see Fig. |J). 

This profile is plotted on a velocity range of -200 to +150 
km s -l < tt Gru C0(2-1) 



bO 





:/s0 



-50 



Fig. 3. RA-declination map of the CO(2-l) emission from 
tt 1 Gru. The offsets are in arcseconds with respect to 
a(1950) = 22 h 19 m 41.13 s , 5(1950) = -46° 12' 02.4". The 
temperature scale for each profile runs from -0.2 to 2 K 
and the velocity scale from -100 to +80 km s . 
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V LSR (km/sec) 

Fig. 4. Dependence of the position of peak CO(2-l) emis- 
sion on velocity (a) upper panel: dependence on declina- 
tion. The ordinate is arcseconds with respect to <5 = — 
45° 11' 00". The error bars are la. The vertical solid lines 
show the velocities of the horns on the CO(2-l) line ob- 
served at the star's position (Fig. |l|). The vertical dotted 
lines show the velocity extent of the SiO(6-5) emission, 
(b) lower panel. As (a), for right ascension. The ordinate 
is with respect to a a = 22 h 19 m 00 s . 



Fig. U reveals no resolvable kinematic structure in the 
EW direction, but shows velocity gradients in the NS di- 
rection, in agreement with S92. Fig. ^ shows that the 
two horns on the CO (2-1) line profile arise from different 
locations in the sky, separated by about 10". 

At higher velocities with respect to the line center, 
±(6 to 12) km s _1 , the declination - velocity curve turns 
over. At higher velocities yet, the position offsets are 
smaller still and decrease essentially to zero. There is a 
second, lower amplitude inflection point in the velocity 
- right ascension curve near the central velocity, with a 
velocity gradient of ~ 5 km s _1 over a region a few arc- 
seconds in diameter. 




18 km s"' 



Fig. 5. Proposed model of the it 1 Gru circumstellar shell. 
The thickness of the model disk is one-fifth of its radius. 
The disk is inclined at 55° to the line of sight. North is up 
in this figure, east-west is perpendicular to the plane of 
the page. This model does not include the fast wind (see 
text) 



4. A Model of the tv 1 Gru Envelope 

First we use a simple model for the molecular envelope 
to estimate an approximate mass loss rate. The SiO(6- 
5) line profile suggests simple spherical outflow, and we 
calculated a model CO emission profile assuming spherical 
outflow at constant velocity and mass loss rate, with a 
relative abundance f = CO/H2 = 6.5 x 10 -4 (Lambert et 

n.-i 



al. |1986|) , V = 10 km s" 1 , and M = 4.5xlO" 7 M Q yr 
The model CO(2-l) line profile is parabolic and does not 
reproduce the observed line shape. 

We then investigated a model consisting of an expand- 
ing flattened system, or disk, tilted to the line of sight. The 
line profiles produced by this configuration were calcu- 
lated from a model of the envelope geometry, with the line 
formation calculated by a Mo nte C arlo radiative transfer 
code based on that of Bernes ( 1979 ) modified for axisym- 
metric geometries. The envelope is modeled by a set of 
concentric rings, with sufficiently fine radial spacing to 
approximate a smooth distribution in optical depth (see 
Crosas & Menten 1997). The rings have 30 radial and 30 
height spacings. Within each of the 30 x 30 rings, the den- 
sity and temperature are constant. The 3D components of 
the velocity vector are calculated at each photon step (cf. 
Bcrncs 1979). The level populations are solved for 30 ro- 
tational levels and for infrared pumping to the first vibra- 
tional state (v = 0— *1 at 4.6 /mi) using the observed 4.6 
/mi flux density of it 1 Gru (Gezari et al. 1993[ ). The mass 
loss rate, radial expansion velocity, turbulent velocity, ki- 
netic temperature distribution and disk radius are input 
parameters. The line profile along a given line of sight is 
calculated by rotating the structure, calculating the line 
emission across the structure and convolving the emission 
with a two-dimensional circular gaussian beam. The fast 
molecular wind was not included in this model. 
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Fig. 6. Comparison of observed and model CO(2-l) line 
profiles for the envelope of tt 1 Gru. The geometry of the 
model is shown in Fig. |^ and described in the text, and 
the calculated line profiles are compared with the data 
for the inner 20" of the envelope (cf. Fig. g). The solid 
line shows the best fit model, which includes a pointing 
offset of +5" in declination. The dashed lines show the 
calculated line profile assuming that the map is centered 
on the center of the envelope (see text). 



A model in which the bulk of the CO(2-l) emission 
arises from a tilted, expanding disk produced by constant 
mass loss fits the data well. This model is sketched in Fig. 
g and the model and observed line profiles are compared 
in Fig. g. The northern part of the disk is tilted away from 
the observer at an angle of 55° to the line of sight. The 
temperature profile is approximated as a power law: T(r) 
= 300 K (r/fO 15 cm) -07 , as found in detailed models of 
AGB s tar w inds (Goldreich & Scoville (|l97§); Kwan & 
Linke ( |l982j ). 

The best fit model has a disk radius of 5 x 10 16 cm, a 
thickness of 10 16 cm, and is produced by a constant mass 
loss rate of 1.2 x 10 -6 M© yr _1 . The radial expansion 
velocity in the plane of the disk is 15 km s _1 , increasing 
slightly to 18 km s _1 at the poles, and the turbulent ve- 
locity is 1 km s~ x . It was also found that better agreement 
with the data is obtained if the telescope pointing is offset 
from the center of the envelope by +5" in declination (note 
that declination and zenith angle offsets are approximately 
equivalent for observations from Hawaii of objects which 
lie at very southern declinations). The outflow velocity in 
the disk agrees well with the value of 13 ± 2 km s _1 ob- 
served for gas in the dense inner regions as observed in 
the SiO(6-5) line (Table 1). 



5. Discussion 

The model shown in Figures 5 and 6, of an expanding disk 
tilted in the north-south direction, reproduces the obser- 
vations well, and also agrees with the higher spatial reso- 
lution observations of S92, which show that the emission 
from the envelope is elliptical with the major axis lying 
east-west. The difference between the discussion by S92 
and that here is one of interpretation; S92 suggests that 
the spatially separated horn features arise from a bipolar 
flow perpendicular to the disk, while our model identifies 
them with the northern and southern halves of the tilted 
disk (cf. Figure 5) whose projected major axis lies east- 
west. While the fast wind from tt 1 Gru is not explicitly 
modeled, the observations suggest that it is likely to be a 
continuation of the velocity increase towards the poles. 

The model disk which reproduces most of the emission 
from the inner envelope is much smaller than the envelope 
extent of~6xl0 17 cm observed via 60 /im emission from 
dust (Young, Phillips and Knapp 1993), and the mass loss 
rate required to produce the observed emission from the 
disk is much higher (1.2 x 10 -6 M© yr _1 ) than the mean 
value (~ 4x 10 -7 M© yr~ x ) given by the CO line intensity 
and 60 fim flux density. The envelope may thus be roughly 
spherical, with a strong density increase towards the plane 
of the disk; alternatively, the mass loss rate of ir 1 Gru may 
have undergone a large increase in the last 1000 years. 

The model proposed here for the envelope of tt 1 Gru 
is similar to that suggested for the envelope of the carbon 



star V Hya by Knapp et al. (1997). The molecular line 
emission from both envelopes is similar: the lines emitted 
from the inner envelope have simple parabolic shapes, the 
CO lines have complex, double horned shapes, and the 
velocity separation of the horns is significantly smaller 
than the velocity range of the "high density" emission. 
Also, both stars have fast (V > 50 km s _1 ) molecular 
winds. Fig. [7] reproduces the IRAS color-color diagram 
for evolved stars with fast molecular winds from Knapp 
et al. (1 997) with the data for tt 1 Gru added Jorissen & 
Knapp ( 1998 ); the colors are plotted for the two epochs 
at which the star was observed. Like V Hya, and unlike 
all the other stars with fast molecular winds, tt 1 Gru still 
has the infrared colors of an AGB star. 

The structure of the tt 1 Gru envelope is remiscent of 
structure seen in some planetary nebulae. For example, 



Manchado et al. (1996) have pointed out the existence of 
a class of quadrupolar planetary nebulae, which they at- 
tribute to the presence of two bipolar outflows ejected in 
different directions. Their proposed model for these flows 



(see also Livio & Pringle 1996; Guerrero and Manchado 
[1998 ) is that the bipolar flow is intermittent and that 
the flow axis precesses due to the precession of a colli- 
mating disk. The similarity of this structure to that sug- 
gested from these observations of the tt 1 Gru envelope 
shows that this complex structure is already present in 
the pre-existing circumstellar envelope. Detailed observa- 
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Fig. 7. IRAS color-color diagram for evolv ed stars with 
fast molecular winds (cf. Knapp et al. 1997 ). The dashed 
lines outline regions I, II, III and VII which are occupied 
by AGB stars. The two connected points for ir 1 Gru corre- 
spond to data taken at different epochs ( Jorissen & Knapp 
1991). 



tions of the molecular emission on angular scales of < 1" 
will become possible before too long; it will be interesting 
to investigate the details of the complex ways in which 
stars end their lives, which are only dimly discernible in 
the present observations. 

6. Conclusions 

1. This paper describes molecular line observations of the 
envelope around the S star ir 1 Gru. We find that the 
envelope contains complex kinematic structure: 

(a) A fast molecular wind, with a maximum outflow 
speed of at least 70 km s^ 1 and perhaps as high as 
90 km s _1 . 

(b) A flattened, tilted disk expanding at 15 - 18 km s _1 
with a mass loss rate of 1.2 x 10~ 6 M Q yr _1 and a 
kinematic lifetime a factor of 10 smaller than that 
found for the circumstellar dust distribution. 

2. Our model of the molecular line emission from the ir 1 
Gru envelope shows that the bulk of the emission arises 
from a disk (and this model also provides a good de- 
scription of the envelope of the carbon star V Hya). 
It is of interest that the scale size of the velocity per- 
turbation in the inner envelope (Figure 4) is similar 
to the projected distance of ir 1 Gru's companion, since 
the separation of the stars is large enough that the 
companion is expected to have very little effect on a 



freely expanding wind (cf. the models of Mastrodemos 
& Morris 1998; see also Jorissen 1998). 
3. The mass loss rate required to produce the flattened 
disk is several times higher than the mean mass loss 
rate inferred from the global CO profile or the 60 /mi 
emission. The mass loss rate of ir 1 Gru may there- 
fore have increased significantly in about the last 1000 
years (the kinematic age of the disk). The star may 
be close to the end of its evolution on the AGB. If so, 
the present observations show that the production of 
a fast molecular wind, and of complex envelope struc- 
ture, begins while the star is still on the AGB. 
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